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anarchic mixing angles are insensitive to the number of singlets, the observed ratios of neutrino 
masses prefer small n for the simplest linear measure. 
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I. INTRODUCTION 

The (type-I) seesaw mechanism [l[ is arguably the simplest and best motivated framework to not only 
give neutrinos mass, but also explain the smallness of said masses compared to the electroweak scale. 
The necessary right-handed neutrino partners are assumed to be heavy to suppress the active neutrino 
masses, which also allows for baryogenesis through leptogenesis via the decay of the heavy states, 
thus explaining the baryon asymmetry of the universe. While initially motivated in the context of grand 
unified theories, the seesaw mechanism has since been applied in its own right. The number of right- 
handed neutrinos, transforming as singlets of the Standard Model (SM) gauge group, is then no longer 
constrained by anomaly cancellations or the need to fill up an irreducible representation, but rather a 
free parameter. At least two right-handed neutrinos are necessary to reproduce the observed neutrino 
mass-squared differences, as well as leptogenesis Q. There is however no upper bound, and models with 
O(10 2 -10 3 ) singlets have been studied in the context of leptogenesis Q, lepton flavor violation (LFV) Q 
and as a way to explain the large mixing angles of the active neutrinos @. Like for most things in 
life, there is also a motivation by string theory Q. The number of singlets is formally infinite in extra- 
dimensional theories, as right-handed neutrinos are not restricted to the SM-brane and will thus lead to 
a tower of Kaluza-Klein excitations H . 

Still lacking, however, is a proper parametrization of the arising mixing and coupling matrices in terms 
of physical quantities, as needed for example for efficient parameter scans (as the number of redundant 
parameters explodes for large n) . The goal of this work is to provide exactly these parametrizations for 
arbitrary n, both for the most general case (top-down approach, to some degree discussed in Refs. @) 
and the seesaw limit (bottom- up). 

The rest of this paper is organized as follows: We first fix our notation in Sec. HH discuss the top- 
down parametrization in Sec. Mil and the bottom- up parametrization in Sec. lIVl also deriving constraints 
on the involved parameters from perturbativity and LFV. To illustrate possible effects of n singlets we 
discuss basis independence (neutrino anarchy) in this framework and present the resulting distributions 
for the neutrino masses in Sec. |Vl App. [X] gives a brief introduction to the group U (N) and various 
representations of its elements, as needed for our parametrizations. App. [B] collects and extends basis 
independent measures for miscellaneous matrices necessary for the discussion of anarchy. Finally, App. [C] 
is devoted to the somewhat special case n = 2, which would interrupt the flow if included in the main 
text, which will only cover n > 3. 
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II. FRAMEWORK AND NOTATION 

Introducing n SM singlet right-handed neutrinos Nj to the Standard Model modifies the Lagrangian 
by the following terms: 



-SM 



iN^Ni - (n 3 {Y v )- a H^L a + ijVj (M R ) jk N£ + h.c.) , (1) 



where a sum over a = e, /i,r (sometimes 1, 2, 3 in the following) and j, k = 1, 2, . . . , n is understood and 
will in the following often be denoted in a vector notation (e.g. as NMrN c ). Without loss of generality 
we work in a basis where the charged lepton mass matrix is diagonal, which can be accomplished by 
unitary transformations of the lepton fields. Mr can be diagonalized by an n x n unitary matrix V 
like V t MrV = diag. This merely redefines Y v , so we will from now on work in a basis where Mr 
is diagonal, with positive real entries M^} After the Higgs doublet H acquires its vacuum expectation 
value v, the Dirac mass matrix for the neutrinos m£> = vY v is generated, leading to a (3 + n) X (3 + n) 
Majorana mass matrix for the neutral fermions (ul,N c ) t : 

Mm = L°_ Tl) • CO 



m D M} 

Assuming Mr 3> mjj gives the low-energy neutrino mass matrix for the flavor neutrinos Vf ~ vl — 

4. -I 

m' D A4 R N c of seesaw form: 

M v ~ -mlM^m D = -v 2 YlM~ R l Y v . (3) 
Diagonalization of _M„ can be performed in the following way 

M v = Up MNS diag(mi, m 2 , m 3 ) ^pmns = ^pmns d ™ ^pmns > ( 4 ) 
with the unitary Pontecorvo-Maki-Nakagawa-Sakata matrix (PMNS matrix) in standard parametrization 

C12C13 S12C13 si 3 e" ,,5N 

f/pMNS = P' ( -C23S12 - S 23 Sl3Cl2e I<5 C23C12 - S 2 3Si 3 Si2e lS S23C13 I P, (5) 

S23S12 - c 2 3Si 3 c 12 e lS -S23C12 - C23Si3«i2e l5 c 23 ci 3 

and the Majorana phase matrix P = diag(l, e m / 2 , e I/3 / 2 ). Here we used the abbreviations Sjj = sin^y and 
Cij = cos Oij for the three mixing angles. The phase matrix P' = diag(e m / 2 , e lb / 2 , e 4C / 2 ) can be absorbed 
by the lepton fields and is hence unphysical. The neutrino mass eigenstates are then v m = Up MNS Vf. 

At low energies the 9 complex entries of the symmetric Majorana matrix Ai v decompose into the three 
eigenvalues m.j (neutrinos masses), three mixing angles #23, #12, ^13 (atmospheric, solar and reactor angle) 
and three CP violating phases 5, a, /3, the latter two being unobservable in neutrino oscillations (but in 
principle testable in Oi//3/3 experiments fToj| ^ . 

III. TOP-DOWN PARAMETRIZATION 

All the information about neutrino mixing is encoded in Y v , the only non-diagonal matrix in the lepton 
sector. For the general complex nx3 matrix Y v , we can write down a singular value decomposition 



/diag(yi,2/2,J/3)\ 
3 



Vl = VrDyVI , (6) 



3 / 



* This is not only convenient, but also necessary to eliminate unphysical parameters. To this effect we note that a diagonal 
Mr can store 0(n) parameters (together with Y v ), while a non-diagonal Mr has 0(n 2 ) independent entries. 
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with positive singular values yt. Out of the 6n real parameters in Y v , 3 phases can be absorbed in the 
lcpton fields, so only 6n — 3 are physical. We will confirm this with an explicit paramctrization of the 
unitary matrices Vr,l below. For the 3x3 matrix Vl, we take the PMNS paramctrization from Eq. ([5]), 
while the n x n matrix Vr can be written as 



n n 



vr = ( n n fi «( a «'fe) I x dia §( £ 

,i=l j=i+l 



e**»). 



(7) 



Here, fiyfoy, £y) denotes matrices of the form 



/I 



V 



cos(a,.,) 



sin(ay)e l ^ iJ ' cos(aij) 



sin(aij)e l ^ i3 ' 



(8) 



V 



where sin^a^e***^ sits at the i-th row and j-th column. See App. [A] for a derivation and a proof for 
the validity of this form of Vr. It is easy to check that Vr (Vl) has the proper n 2 (3 2 ) real parameters, 
n(n + l)/2 (6) of which are phases. From the form of Y v in Eq. ©, it is clear that the angles <j>j either 
act on zeros in Dy or can be absorbed by the phase matrix p£ in v[. An overall phase and P[} can 
be absorbed by the lepton fields, so Vl contains only 3 angles and 3 phases. The ordering of the tlij in 
Vr is of course arbitrary, which allows us to move all the rotations that act on the zeros in Dy to the 
right, so (n — 3)(n — 4)/2 matrices of the type (O drop out, leaving n(n — 1) — (n — 3)(n — 4) = Qn — 12 
real parameters in Vr, half of which are phases. Consequently we can restrict the product in Eq. |(7J| to 
n»=i rij=i+i- § A s expected, Y v contains 3 (from yi) plus 6 (Vl) plus 6n — 12 (Vr) equals 6n — 3 real, 
physical parameters, 3(n — 1) of which are phases. This agrees with the parameter counting in Ref. pj| . 
The full mass matrix can then be written as 



■Mfull = 







WtDTyVl 



vVrDyVI Mi 



(9) 



with the explicit form for Vl'- 



V l 



L „L 



CinC 



L „L 



C90S 



23 s 12 



12^13 

„i „L „L ip£ 
i 23 A 13 c 12 e 



S 23 S 12 c 23 S 13 C 12 eJ 1 



r L r L - 
C 23 L 12 

6 23 c 12 



„L L 
6 12 c 13 

_ „£ „L L -ifi^ 



, 23 ; '13' : '12 ,: 
c 23 s 13 s 12 e 



i 23 c 13 
C 23 L 13 



diag(l,e^/ 2 ,e^/2 ) . 



(10) 



and the shorthand notation sf. • = sin fl^- etc. for the mixing angles. In the seesaw limit, we thus find 



M v si -v 2 VlD^(V ] R M R V^)- l D Y Vl 



(11) 



Note that the occurring mass matrix V^MrV^ is for n > 3 no longer the most general complex symmetric 
n x n matrix, because we already restricted Vr to the physical angles. 

From Eq. ([9]) or Eq. (Ilip one can now extract the physical masses and mixing angles in terms of the 
parameters yi , ay , £y , 8fj and /?f . 

For completeness, let us count the total number of parameters in the lepton sector. We have 3 charged 
lepton masses, n heavy neutrino masses and 6n — 3 real parameters in Y u , summing up to In real 
parameters. 



§ In other words, the U(n — 3) subgroup of U(n), acting on the lower right parts, can be modded out. 
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IV. BOTTOM-UP PARAMETRIZATION 



Seeing as the seesaw formula (jlll) contains more parameters than measureable observables, it proves 
convenient to parameterize our ignorance of the high-energy sector of the seesaw model. In other words, 
there are infinitely many different Y v and M.r that lead to the same low-energy neutrino parameters 
via the seesaw formula M. v ~ —rnj^JA^mD- The parametrization of this ambiguity is the goal of the 
Casas-Ibarra-parametrization (l2j of the Yukawa coupling: 



Y v = - y /M^Ry/d^Ul lSNS . (12) 



,, . . A ™ ^pmns ■ 

This Y v solves 

— v 2 Y^ M j^Y '„ ~ M u = ?7p MNg d m Up MNS , (13) 

as long as the complex n x 3 matrix R satisfies R T R = tsx3- Correspondingly, any values of Mr and 
R in Eq. (1121) lead to the same low-energy model, which means we can fix the PMNS mixing angles and 
masses to the observed ones and study Y v independently in other processes. 

We will now give an explicit parametrization of i?, which has 3n complex entries obeying 12 real 
constraints (R T R = 1 is symmetric); we therefore expect to find 6(n — 2) real physical parameters. This 
agrees of course with the parameter counting in the last section, as we also have 3 charged lepton masses, 
3 light neutrino masses, n heavy neutrino masses and 6 parameters in the PMNS matrix, summing up to 
7n input parameters. Since we only need three orthonormal complex vectors £ C™ to define R, we can 
write 



R = O n 



/1 3 X 3 \ 

3 



VOW 



C„xnlnx3, (14) 



with a complex orthogonal n x n matrix O nxn £ 0(n, C). Similar to the U(n) parametrization used for 
Vr in the previous section, we can write O nxn as a product of the (7^ , £L ) matrices from Eq. |8]). 
but with the arguments £L = and 7^ £ C in order for O nxn to be complex orthogonal: 

n n 

O n xn = S X 

n n n «(7«,o). (15) 

i=l j=i+l 

Here, S is a diagonal matrix with entries ±1, which decides the determinant of O n xn and therefore 
whether it describes a rotation or a reflection. We can once again mod out the 0(n — 3,C) subgroup 
that acts on the zeros in l nX 3 by reordering the rotations, so the product in Eq. (|15p can be restricted to 
Yli—i rifci+i- Generally speaking, R is an element of the quotient space 0(n,C)/0(n — 3, C). For real 
angles 7y, R is an element of the compact Stiefel manifold V^IR") = 0(n,M)/0(n — 3,R), which admits 
a Haar measure and could thus be used for statistical considerations. 

Note that the case n = 4 is special in the sense that there are no superfluous rotations that need to 
be modded out. One can therefore take any complete parametrization for O^xi without having to worry 
about unphysical parameters. 

For complex angles, O nX n and therefore R and Y v can have arbitrarily large elements, because the 
entries cos ix = cosh a: grow exponentially with a;. As a result, the Yukawa couplings Y „ can be very large 
even for small A4r, which is somewhat counterintuitive to the seesaw formula and relies on cancellations 
in the matrix product Y v M^Y v . Even though these cancellations suggest finetuning of parameters, this 
is a viable possibility and deserves discussion, especially because a low M.r with large Yukawa couplings 
makes an LHC discovery (direct or indirect) possible [13j. Popular models to obtain this finetuning in a 



natural way are for example inverse- [14| and linear seesaw 15j , usually based on the case n — 6 with an 
imposed structure in M. r. and Y v by some symmetry. 

A useful constraint on the complex mixing angles 7^ comes from the perturbativity of Y '„, on which 
many calculations involving Y v rely on. For this, we propose yet another, slightly different parametriza- 
tion of -R, which we obtain by noting that f2y (7^ , 0) = fiy (Re (7^ ) , 0) x J7y (ilm (7^), 0), and that the 
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ordering of the f2,j is arbitrary: 

r=sx (n f[ Qijivij^u (n n m^o)]!^. m 

\i=l J"='t+1 / \i=lj=i+l / 

The 6(n — 2) parameters that make up R are thus equally divided among the new real parameters rjij and 
Pij. The product involving rjij is just an orthogonal matrix, so we can restrict rjij e [0, 2tt) or [— 7r, tt). " 

The parametrization (fTfJ|) can be more useful than ([Til) , (fT5|) because the 77^ drop out of the expression 
R)R, which arises from (the often occurring) Y*JY V for degenerate Mr. 

A. Perturbativity 

Perturbativity of the Yukawa couplings gives constraints of the form [ItJ 

3 



1 

tr(Ft^) = -J2J2\^\ 2M ^ % °(!) ■ 

This inequality is satished if 



i=i j=i 



3 /ITeV lO- 2 eV 



^^Tinjsr^ 10 ^^— V— • (18) 

which constrains each entry in _R. Note that the upper bound goes down for increasing n. Once the 
smallest masses are specified, the above equation can be used to give a total upper bound on the param- 
eters pij. For example, the largest entry in R can be estimated as coshp max if just one is large, and 
(cosh/9 max ) 3 (™~ 2 ) if all p^ are large. This typically leads to upper bounds on \pij \ of 0(1-10). 



B. Lepton Flavor Violation 



In the seesaw limit, the neutrino states vr, that couple to the charged current are approximately given 
by [3 vl ~ (1 — ^m^M^mD^PMNS^m + mtM^ 1 N c , with the mass eigenstates v m and iV c . This 
leads to lepton flavor violating decays I — > £'j [19|, of which p, — »■ e 7 gives the strongest constraint: 



BR(^i e 7 ) 



3a 
8^ 



R 2 m,D C^pmns B C/ PMNS + h.< 



(19) 



with the diagonal matrices Aij = Sij g(Mj /M^) and Bij = 5ij g(m 2 j / M^) and the monotonic decreasing 
loop function 

10 - 43a; + 78a; 2 - 49a; 3 + 4.x 4 + 18a; 3 log a; 



9(x) = 



6{x- l) 4 



9(0) = 



5(1) 



17 
12 



j(x-kx>) = -. (20) 



Here wc already omitted the contribution from the light neutrinos, as the corresponding term 
(C^pmns B U PMNS )fj, e is suppressed by the tiny mass-squared differences of the active neutrinos The 
remaining contributions are naively of order (mu/yVfi?) 2 , but can be enhanced via R. In the approxima- 
tion Mj 3> My/ ^> mfc we have A oc 1 oc B, so Eq. (|19l) simplifies to 



BR(p -4- e 7 ) 



3a 
8^ 



3 a: 
8^ 



(u PM mV^R t M B 1 R\/d^Ul 



PMNS 



(21) 



' The parameter space of the angles can be further reduced using discrete field redefinitions like in Refs. Il6l . which however 
lies outside the scope of this letter. 
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Mi/Mo M t /M„ 

Figure 1: Distribution of the Mr eigenvalues Mk according to the linear measure from Eq. ()23|) with the boundary 
constraint YJl=i M l ^ M o ■ 

Since this branching ratio can easily exceed the current limit BR(/i — > ej) < 2.4 x 10 -12 [2(| — for 
perturbative couplings — it gives constraints on the pij. Naively, the relevant product in Eq. (j2"T]) is 
\R\ 2 m/M, so the LFV constraints have a different structure than the perturbativity bound (JT5J) . Note 
that this constraint is connected to the unitarity of the lepton mixing matrix [2l|, because the relevant 
matrix (1 — ^rn i D A4] i 2 rnD)UpMNS is no longer unitary at this order. Consequently, the above bound 
also ensures the validity of our parametrization for Y v from Eq. (|12p. as the matrix C/pmns m Y v only 
corresponds to the lepton mixing matrix if (m D Ai] i 2 rnD)ij *C 1. 

While not particularly nice to look at analytically, the constraints given here should be useful in 
numerical scans involving Y v for arbitrary many right-handed neutrinos, extending for example the 
analysis of Ref . [TtJ ■ 



V. APPLICATION TO NEUTRINO ANARCHY 



The idea that the large lepton mixing angles and small neutrino hierarchy are not due to some flavor 
symmetry but rather the absence of any distinction between the neutrino generations has been around 
for over a decade [22], |23| . The basis independence of the neutrino mass and mixing matrices leads to a 
distribution of the mixing angles according to the Haar measure of the Lie group G that diagonalizes M v 
(so G — U(3) for complex, G — 0(3) for real M v ). Since these distributions prefer large mixing angles, 
this ansatz seems to fit well to the observed large leptonic mixing angles (see however also a critique of 
the anarchy approach in Refs. [24], Hfl). While the distribution of the mixing angles is uniquely given by 
the Haar measure, the distributions of the masses and Yukawa couplings are not unique, and usually the 
simplest (linear) measure is used (listed in App. IB]) . 

Let us now discuss neutrino anarchy, i.e. basis independence, with n right-handed neutrinos. In the 
seesaw limit we have the low-energy mass matrix for the active neutrinos from Eq. 

M v * -v 2 VID^M^VrDyVI , (22) 

where Vl and Vr are now assumed to be distributed according to the Haar measure of U(3) and U(n), 
respectively, while we take the simplest linear measure for the eigenvalues of Mr [23| 

n n 

dM R cc 11 \Mf - Mf\ 11 M k dM fc . (23) 

i<j k=l 

We employ the invariant boundary tr(M R MR,) = < Mq for the scanning region; see Fig. [T]for a 

visualization of the distribution.** For the singular values of Dy we use the linear measure (see App. iBl 



** On a technical note, random unitary matrices following the Haar distribution were obtained from a QR decomposition |2(| . 
while random entries out of the linear measures dMn and dDy were drawn using a Metropolis algorithm |27| . This 
algorithm is fast even for large n and conveniently works without knowing the normalization factor of the distribution. 
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for the derivation) 

3 3 3 

dD Y cx J] (yl) n - 3 f[(yf - y]f f[ y k dy k , (24) 

m=l i<j fe=l 

with a similar boundary tx{Y\ / Y v ) = ^(DyDy) = Dj ^ Vo- The effect of large n will be a reduced 
hierarchy in the singular values yi (see Fig. [5]). Note that Eq. (pM)) is only valid for n > 3, see App. [Clfor 
the used measure for n = 2. 

The value v 2 y$/MQ will fix the overall light neutrino mass scale and can be used to fix one of the 
mass-squared differences Am|. The PMNS mixing angles and phases do not change with n and are 
given by the Haar measure for 17(3) [23| (ignoring unphysical phases): 

dC/pMNS oc ds^ 2 ds^ d c i3 d<5 da d/3 , (25) 

so the only effect of n ^ 3 is a change of the three eigenvalues of Ai v , i.e. the ratio R v = (m| — rri\)/ {m\ — 
mf), where we sorted the masses like m\ < iri2 < m^. In this notation, R v < 1/2 corresponds to normal 
hierarchy (NH) and R v > 1/2 to inverted hierarchy (IH), with best-fit values R^ H = 0.03 and R™ = 0.97 
(IH), respectively, taken from a recent global fit [28j . 

We observe from Fig. [3] that for increasing n, the distribution for i?„ shifts to larger values, while the 
width decreases. This reduction of hierarchy can be tracked back mainly to the Yukawa couplings, as the 
factor nL=i(2/m)"~ 3 m d-Dy pulls the yi tightly together. An analysis without this factor shows that 
the behavior of R v goes in the opposite direction, i.e. the maximum shifts towards small values for large 
n while being diluted. The distribution then however quickly converges as the only n dependend change 
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comes from VrMrVr, for which only the first dozen or so Mj are relevant, due to the hierarchy in 
Mr. Correspondingly, to shift the maximum of i?„ to its observed value, one not only has to omit the 
factor Y\ m=1 (y m ) n ' 3 m ^Dy, but also increase the hierarchy by inserting something like Y\i<j{vl ~ v]) 2 ' 
i.e. using a non-linear measure (see Ref. j| for other yi distributions that accomplish this task). We 
conclude that the large n limit in the anarchy approach worsens the agreement with data. While the one 
data point Nature provides can obviously not be used to find the overlying distribution, it should be fair 
to say that anarchy works best for small n. 

Let us finally comment on the limit n —> oo. As the Yukawa couplings yi become quasi-degenerate 
for very large n, the neutrino masses are simply given by the upper left 3x3 submatrix of V^Ai^ Vr 
(times a prefactor). Since there is no preferred direction in the Vr rotations, the entries (Vr)^ have the 
same mean magnitude for a given n. The mean entries ((VjMj vR)jk) then take the form of a random 
walk with decreasing step-size J2 m e%dm /M m , and due to the hierarchy in the singlet masses Mk, only 
the first couple of steps are relevant. The mean magnitudes are again the same for all entries, so the only 
structure in the neutrino mass matrix comes from random phases Ojk'- ((Vr M^Vr) jk) oc e t6jk . It is 
easy to show that a matrix of this type predominantly yields NH solutions, so even though IH becomes 
more and more probable for large n, it will never dominate. 

VI. CONCLUSION 

One of the simplest ways to explain the masses of the active neutrinos is the introduction of right- 
handed partners. As the number n of these SM singlets is in principle unconstrained, it is interesting 
to study implications of varying n. We have given explicit parametrizations for the mixing and coupling 
matrices connecting the n right-handed neutrinos to the SM in terms of physical parameters, both in 
the top-down and the bottom-up approach (Casas-Ibarra-parametrization). For the latter, constraints 
on the involved parameters from perturbativity of the Yukawa couplings as well as charged lepton flavor 
violation have been discussed. 

As a novel application of the n singlet framework we studied basis independence in the neutrino sector, 
i.e. anarchy. Of the low-energy neutrino parameters, only the neutrino mass distribution changes with n, 
and we showed that anarchy — with the simplest linear measure — seems to prefer small n in view of the 
observed mass-squared differences. 
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Appendix A: The Group U(N) 

For convenience we provide a short and rather colloquial review of the group U(N, C) = U (N) . Formally 
one can define it as the set of complex N x N matrices that leave the inner product of C N invariant, 
i.e. U(N) = {V € C NxN | VW = 1}. With a little effort one can further show that the entries of U(N) 
form a group under matrix multiplication, actually even a compact connected Lie group. Counting the 
number of entries in V £ U(N) and constraints, we find that U(N) has N 2 real parameters. One can 
show that every element in U(N) can be written as V = cxp(A) with skew-hermitian A = — A* . We can 
choose a basis for A = such that the N 2 linearly independent Xj satisfy Lie algebra relations 

[X k ,Xi] = ifkimXm with structure constants f k i m - 

The N 2 generators Xj can be further separated into N(N + l)/2 symmetric generators X s — iY — iY T 
and N(N — l)/2 antisymmetric real generators X a = — Xj. If also proves convenient to treat the N 
diagonal (symmetric) generators Yd separately, as they are especially easy to exponentiate. To fix the 
normalization, we chose (Y^)^ = SijSjk, while the N(N — l)/2 non-diagonal Y s and X a have just two 
non- vanishing entries, namely +1 in the upper right half and ±1 in the lower left (plus sign for Y s , 
minus sign for X a ). Note that this normalization — and therefore also fkim — differs from the common 
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convention, but is of no importance in the following (see also Refs. [9j). A general element of U(N) can 
then be written as the product of all rotations: 



N 



N(N-l)/2 N(N-l)/2 

U(N)3V = HeM^ 3 Yj) Y[ eM^kY s k ) J] exp( 7m X") . 



(Al) 



j'=i 



fe=i 



The equivalence of this parametrization to the initial V — exp(A) can be proven using the Baker- 
Campbell-Hausdorff formula and the fact that the generators satisfy a Lie algebra. The order of the 
rotations in Eq. (jAll) is of course arbitrary [29[ and should be chosen to simplify given expressions. The 
three different types of rotations take the form 



A 



\ 



cos(/3) isin(/3) 



i sin(/3) cos(/3) 



/I 



cos(7) sin(7) 



sin(7) cos(7) 



and 



V 



exp (iaYd) — diag(l, . . . , 1, exp(ia), 1, . . . , 1) , 



(A2) 



V 



(A3) 



with real angles a, j3, 7. Another useful parametrization for a general unitary matrix along the same lines 
as above is given by 



N 



N(N-l)/2 

U(N) 3V = H exp^.Fj) J] eM^' k Y s k + i k X k a ) . 
j=l k=i 



(A4) 



where we combined the generators X a and Y s that have the same vanishing entries. This yields a product 
of matrices of the type 



A 



\ 



cos(zj) sm(zj)e l £ j 
■ sm(zj)e~ l tj cos(zj) 



V 



(A5) 



1/ 



with the real parameters Zj — \v,j\, ^ = arg(uj) and Uj = i/3'j — 7^ in our chosen normalization. Once 
again, the ordering of the matrices in Eq. (|A4|> is arbitrary. 



Appendix B: Measures for Anarchy 

This appendix provides the necessary measures for the anarchy approach in Sec.[V] We refer to Ref. [23| 
for a detailed derivation of the known measures, and simply quote the result for real and complex 
Majorana N x N matrices: 

N N N N 

dMcai oc JJ I A - A I II dDkd ° ' decomplex cx JJ |U? - £>f I D k dD h dU , (Bl) 

i<j fc= 1 i<j k—1 
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where the decompositions M Ica .\ — ODO T and A'lcompiex = UDU T were used. dO and dU denote 
the usual Haar measure for O(N) and U(N), respectively, while the diagonal matrix D contains the 
eigenvalues of Ai. We will omit the calculation of the normalization factor for all measures in this 
appendix, as they are irrelevant for our considerations. It can in principle be calculated by integration 
over some volume, e.g. tr ( M. l nm plfiY M. complex) < Mq, up to the arbitrary overall mass scale M . 
For Dirac matrices, the measures read [23[ 

N N N N 

dm rcal cx J] \D- - D) | Y[ dD k dO L dO R , dm complcx cx J](£>? - D]f J] D k dD k dU L dU R , (B2) 

i<j k—1 i<j k—1 

with m roa i = OrDO\, m comp i ox = UrDU^, and the Haar measures for O(N) dC>L,R and U(N) dUL,R. 
respectively. For dm com piex, unphysical phases in the overlap of dUi, and dUn should be modded out, 
once a parametrization is specified. 

These invariant measures are of course not unique, but rather the simplest consistent ansatz that give 
the right mass dimension.^ With the above equations, one can easily write down the measure for an 
hermitian matrix H = UhW = W under U(N): 

N N 

dH cx Y[(hi - hjf Y[ dh k dU . (B3) 

i<j fe=l 



As a cross check, one can plug the hermitian matrix "T-compicx mcom pi cx ^q- (jB3|) to rederive Eq. (IB2|) . 
We can now give a measure for a non-quadratic n x 3 complex Dirac matrix M$vd — UrDyU^ with 
Ur £ U(n), Ul e U(3). Considering mI vd Msvd shows that dMsvD should at least contain a factor 

3 3 

H(yl - v)f II y^yk , (B4) 

i<j k=l 

which however has too small of a mass dimension. The same procedure for the n x n hermitian 
MsvdM1 vd runs into problems, because n — 3 of the eigenvalues are zero and would give a vanish- 
ing measure. Taking the product only over the non-vanishing factors however yields 

3 3 3 

dMsvD « II (y-)"~ 3 II(^ - II Vkdy k dU L dU R , (B5) 

m=l i<j fe=l 

which has the right mass dimension (6n) and can therefore be viewed as the linear measure for M$vd- The 
overlap of dUh and dt/j? — this time including angles, see Sec. lIIU — should again be modded out in a given 
parametrization to avoid overcounting of physically equivalent configurations (similar to gauge-fixing, as 
pointed out in Ref. |25[). 

The measure (|B5[) can of course be easily extended to a general measure over 0{N) x 0{M) or U(N) x 
U(M) with M < N: 

M M M M 

dM rcal cx [] \y1 - y) | J] V l - M dy k dO L dO R , dM complex cx JJfc/? - y]f J[ yl +2(N ~ M] >dy k dU L dU R , 

i<j k—1 i<j k—1 

(B6) 

with the N x M matrices M rea i = OrDyO\ and Af comp i ox = UrDyUI. Unphysical rotations need to 
be modded out again. A simple check shows that the mass dimensions are right and N = M indeed 
leads back to Eq. (|B2|) . Eq. (IB6|) can actually be obtained from exactly these requirements, without 



tt For example, -M comp iex contains N(N + 1) independent real parameters, so the measure dMij has mass dimension 
N(N + 1). Rewriting dAiij in the form (IB 1 1) matches this mass dimension, as can be easily verified. 
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the need to use hermitian matrices in the derivation. For example, dMcompiex should contain the known 
N = M factor Ylu:j{Vi ~ vf) 2 Tlk=i ^Vk' which leaves some factor of mass dimension 2M to the power 
of N — M. Since a factor Y[i<j(yi ~ Vj) 771 — which would increase the hierarchy in the yt — has mass 
dimension M(M — l)km/2, it can only be used in special cases, not for general M (at least for integer 
k and m). Thus, for general M, the only factor with mass dimension 2M is Yii ufi leading back to 
Eq. l|B6p. 



Appendix C: The Case n — 2 

In this appendix we collect the formulae for the case n = 2 (see also Refs. 

dUl), which is a little 

different than the n > 3 cases and would interrupt the flow if included in the main text. With the 2x2 
diagonal matrix Mr we find the top-down parametrization in the singular value decomposition for the 
Yukawa coupling 

^=Mo\°£K < ci > 

Two phases of Vr can be absorbed by v£, three of Vl in the lepton fields. Since one of the diagonal 
phase rotations of Vl acts on the zeros, there are only two phases in Vl left, so we can write 

/ cos(ai 2 ) sin(ai 2 )e^ 12 \ ,„ . 



Vl 



— sin(ai2)e '^ 12 cos(«i2) 

L 12 c 13 A 12 L 13 *13 e 



-C 



L „L L L L i&L L L _ L L L iff L L diafffl e^ 2 /2 1) (C3) 

23 s 12 s 23 s 13 c 12 e c 23 c 12 s 23 s 13 s 12 e s 23 c 13 I uld, 6l 1 ; e i L ) > \^°) 



\ S 23 S 12 c 23 S 13 C 12 el ' 31 S 23 C 12 ^^l^Vl^ 1 C 23 C 13 

with the shorthand notation — smBfj etc. for the mixing angles. The linear measure for 2/1.2 for the 
anarchy framework can be derived along the same lines as in App. [B] with the result 

AD Y cx (yl - V 2 2 ) 2 vlvlA Vl Ay 2 . (C4) 

Finally, the bottom-up parametrization is give n by Eq. (|12D with two different R matrices depending on 
the low-energy neutrino hierarchy (7 G C) 



givei 

m 



R^=(° n C ° S7 " Sin7 ) , Rm = f C ° S7 " Shl7 ■ (C5) 

VO sin 7 COS7 I Vsm7 COS7 01 y ' 
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